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Using about 23M BB events collected in 1999-2000 with the BABAR detector, we report the obser- 
vation of several hundred B — > D^*'D K decays with two completely reconstructed D mesons. 
The preliminary branching fractions of the low background decay modes B° — > D*~ D^*'° K + are 
determined to be B(B° -» D*-D°K + ) = (2.8 ± 0.7 ± 0.5) x 10~ 3 and B(B° -> D*-D*°K + ) = 
(6.8 ± 1.7 ± 1.7) x 10~ 3 . Observation of a significant number of candidates in the color-suppressed 
decay mode B + — > D* + D*~ K + is reported with a preliminary branching fraction B{B + — > 

43 ■ D* + D*-R' + ) = (3.4 ± 1.6 ± 0.9) x 10~ 3 . 
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1 Introduction 

Decays of B mesons that include a charmed and an anti-charmed meson are expected to occur 
through the b to c quark transitions b — > cW + , where the 1U + materializes as a cs pair. These 
transitions are responsible for most of the D s production in B decays. D s production has been 
thoroughly studied in experiments running at the Y(4S) resonance [jl], ^, ||]. The inclusive rate 
for D s production in B decays was recently measured by BABAR, where a preliminary branching 
fraction pj: 

B(B - D S X) = (10.93 ± 0.19 sta4 ± 0.58 SJ/st ± 2.73^)% 

is reported. 

Until 1994, it was believed that the cs quarks would hadronize dominantly as D s mesons. 
Therefore, the branching fraction b — > tics was computed from the inclusive B — > D s X, B — > 
(tic) X and B — » S C X branching fractions, leading to B(b — > ccs) = 15.8 ±2.8% ||. Theoretical 
calculations are unable to simultaneously describe this low branching fraction and the semileptonic 
branching fraction of the B meson ||. It has been conjectured that B(b — > tics) is in fact 
larger and that decays B — > DDK (X) (where D can be either a D° or a D ) could contribute 

significantly. This might also include possible decays to orbitally-excited D s mesons, B — > D D**, 
followed by £>** -» £>M K. 

Some experimental support for this picture has been published in the last few years. The most 
significant results are the evidence for wrong-sign D production in B decays (CLEO), yielding 
B(B — > D X) = 7.9 ± 2.2% []|], and the observation of a small number of completely reconstructed 
B -» D^D^K decays, by both CLEO § and ALEPH 0. 

1? — > D^*'D^*'K decays can occur through three different processes: pure external diagrams, 
pure internal (color-suppressed) diagrams and the sum of both. Fig. |l] shows the three possible 
types of decays for charged and neutral -B's. 

In BABAR, the high statistics available allow comprehensive investigations to be made of the 
b — ► tics transitions. In the analysis described in this paper, we present evidence for the decays B — > 
£)(*>D Xg and B — > D^*'D K , using events in which both D's are completely reconstructed. 
After describing the data sample and the event selection, we show the D^*'D K signals for the 
sum of all B submodes. The branching fractions for some of the cleanest modes, such as B° — > 
D*~D^*' K + , are computed and the main systematic errors are discussed. Observation of several 
candidates in the color-suppressed decay mode B + — > D*~ D* + K + is also reported. 

2 The BABAR detector and dataset 

The study reported here uses 20.7 fb~ of data collected at the Y(4S) resonance with the BABAR 
detector, corresponding to (22.7 ± 0.4) x 10 6 BB pairs. 

The BABAR detector is a large-acceptance solenoidal spectrometer (1.5 T) described in detail 
elsewhere |ll]]. The analysis described below makes use of charged track and 7r° reconstruction 
and charged particle identification. Charged particle trajectories are measured by a 5 layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH), which also provide ionisation 
measurements (dE/dx) used for particle identification. Photons and electrons are measured in the 
electromagnetic calorimeter (EMC), made of 6580 thallium-doped Csl crystals constructed in a 
non-projective barrel and forward endcap geometry. Charged K/ir separation up to 4GeV/c in 
momentum is provided by a detector of internally reflected Cherenkov light (DIRC), consisting of 



12 sectors of quartz bars that carry the Cherenkov light to an expansion volume filled with water 
and equipped with 10752 photomultipliers. 

3 Analysis strategy 

The B° and B + mesons[] are reconstructed in a sample of multihadron events for all possible DDK 
modes, namely B° -► D^~D^ K+, D^'D^+K , D { * )0 D^°K° and B+ -» D ( * )0 D^+K°, 
D { * )0 D^ K + , D { * )+ D^-K+. 

The Kg candidates are reconstructed from two oppositely charged tracks coming from a common 
vertex displaced from the interaction point by at least 0.2 cm and having an invariant mass within 
±9 MeV/c 2 of the nominal K° mass. The 7r° candidates are reconstructed from pairs of photons, each 
with an energy greater than 30MeV, which are required to have a mass 115 < M 77 < 150MeV/c 2 . 
The 7T° from D*° must have a momentum 70 < p*(77) < 450MeV/c in the T(4S) frame, while the 
7T° from D° — ► K~7t + -k° must have an energy E(7T°) > 200 MeV. Finally, a mass-constraint fit is 
applied to all the 7r° candidates to improve the energy resolution. 

The D° and D + mesons are reconstructed in the modes D° — ► K~ir + , K~ir + ir°, K~Tr + ir~-K + 
and D + — > K~ir + ir + , by selecting track combinations within ±2<r or ±3(7 of the nominal D mass, 
where a is the mass resolution for the D decay mode considered and the tighter la mass interval is 
applied for B modes with a larger combinatorial background. The K and ix tracks are required to 
be well reconstructed in the tracking detectors and to originate from a common vertex. Charged 
kaon identification, with information from the Cherenkov angle in the DIRC and from dE/dx 
measurements in the drift chamber and in the vertex detector, is required for most D decay modes, 
as well as for the K from 5's. 

D* candidates are reconstructed in the modes D* + — ► D°tt + , D*° — ► D°tt° and D*° — ► D°7, 
by combining a D° candidate with a 7r~, 7r°, or photon. A ±3er interval around the nominal 
AM = M(D*) — M(D°) mass difference is used to select D*'s. Partial reconstruction of D*°'s (no 
7T° or 7 reconstruction) is also used in the B° — > D*~ D*° K + mode, as explained below. 

B candidates are reconstructed from the D^*\ D and K candidates. A mass constraint is 
applied to all the intermediate particles (D°, D~, Kg). Since the B mesons are produced via e + e~ 
—> T(4S) — > BB, the energy of the B in the T(AS) frame is given by the beam energy E£ eam , which 
is measured much more precisely than the energy of the B candidate. Therefore, to isolate the B 
meson signal, we use two kinematic variables: AE, the difference between the reconstructed energy 
of the B candidate and the beam energy in the center of mass frame, and mES) the beam energy 
substituted mass, defined as 



m ES = V^eam " P*B 

where p* B is the momentum of the reconstructed B in the T(AS) frame. Signal events will have AE 
close to and mES close to the B meson mass, 5.729 GeV/c 2 . When several candidates are selected 
per event in a specific B submode (e.g. B + — > D°D K + ), a x 2 value, taking into account the 
difference between the measured and the PDG values of the D masses and of the AM (for D*'s) 
is constructed and only the candidate with the lowest x 2 value is kept for the given submode. 



1 Charge-conjugate reactions are implied throughout this note. 



4 Evidence for signal in the sum of all B submodes 

We present here the distributions obtained by summing all possible B — > D^*'D^*'K decay channels, 
for neutral and charged B decays respectively. 

Since multiple candidates are removed only submode by submode the same event can appear 
several times in distributions obtained by summing over all modes. In the AE distribution this 
manifests itself as three distinct peaks, as can be seen in Figs. || and y. An event will appear in 
the peak near OMeV when reconstructed correctly, in the peak around — 160 MeV if it is a D*DK 
(D*D*K) decay reconstructed as DDK (D*DK), and near the peak around +160 MeV if it is a 
DDK (D*DK) decay reconstructed as D*DK (D*D*K). 

About 120 B 0, s and 180 B decays have been reconstructed. The m^s distributions (Figs. |2| 
and ||) contain only events with \AE\ < 24 MeV. From Monte Carlo studies, the tties resolutions 
of the different sub-modes are quite similar and the 771es spectrum of B° and B events can be 
fitted by the sum of a background shape and a Gaussian function used to extract the number of 
signal events. The background is empirically described by the function 




dN I mL 
oc m ES X W 1 - f J x exp 

dm ES V E beam 

where the only free parameters are £ and the normalization factor. This function is referred to as 
the ARGUS function in the following. The AE distributions (Figs. [2| and ||) contain only events 
with mES > 5.27 MeV/ c 2 . They have been fitted by the sum of a polynomial background and three 
Gaussian functions for the three signal components described above. However, the fits to the AE 
distributions are only indicative since they merge many B and D sub-decay modes, which have 
significantly different AE resolutions depending on the number of 7r°'s or photons involved 

5 Measurement of exclusive branching fractions 

In this section, we present measurements of branching fraction for the three decay channels B° — * 
D*-D°K+, B° -» D*-D*°K + and B + -» D*-D* + K + . Several candidates are also observed in 
the CP conjugate modes B° — ► D^~ D^*' + Kg but without extracting branching fractions. 

5.1 Monte Carlo samples and efficiencies 

The selection efficiencies for each mode were obtained from detailed Monte Carlo simulation, in 
which the detector response is modeled with the GEANT3 program |0|. In addition, data was used 
whenever possible to determine detector performance and the simulation adjusted accordingly. B 
meson decays to DDK were generated with a three-body phase space model. For each sub-decay 
mode, samples of 5000 signal events were produced. Typical efficiencies range from 10%, for 
B° -> D*-D°K + with both L>°'s decaying to Kit, to less than 1%, for B + -» D*-D* + K + with 
D s decaying to Ktttt or K3n. 

5.2 Systematic uncertainties 

Systematic errors account for the uncertainties on tracking and ir° reconstruction efficiencies, K 
identification efficiency, D and B vertexing requirements, efficiency of the requirement on AE used 
to define the signal box, efficiency of the D mass requirement; uncertainty on the background 
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shape; uncertainties on the D and D* branching fractions; uncertainties on the selection efficiencies 
arising from Monte Carlo statistics; and uncertainty on the number of produced BB events in the 
data sample. The breakdown of the different contributions to the systematic error for each mode 
is given in Table |l|. 

5.3 B + -> D*-D* + K + 

The toes distribution obtained for events with \AE\ < 24MeV is shown in Fig. || for the sum of 
all the six possible D° x D decay combinations. A fit to the data is performed with the sum of 
a Gaussian function for the signal and an ARGUS function for the background. The number of 
signal events is 8.2 ± 3.5 and the number of background events given by the ARGUS function is 
1.7. The probability that the signal arises from a background fluctuation is 1.4 x 10~ 5 (> 5a). The 
corresponding preliminary branching fraction is measured to be 

B{B + -► D*-D*+K + ) = (3.4 ± 1.6 ± 0.9) x 1(T 3 

The first error quoted is statistical and the second is systematic. The different contributions to the 
systematic error are given in Table |l|. 

5.4 B° -► D*~D^°K + 

In this analysis we require that either the D° or the D decays to Kir and we do not explicitly 
reconstruct the tt° or the photon from D*° — > D°ir or D 7. The toes versus AE distribution of 
D*~ D° K + combinations is shown in Fig. |B| for the sum of the three D°D sub-modes considered. 
Despite the background level, two separate enhancements, due to the decay modes B° — > D*~ D°K + 
and B° — > D*~ D*° K + , are clearly visible. The enhancement in the region AE ~ 0, toes — 
5.28GeV/c 2 corresponds to decays B° — ► D*~D°K + , while the second enhancement in the region 
AE ~ -154 MeV, toes - 5.28GeV/c 2 corresponds to decays B° -► D*-D*°K+. 

Events containing B° — > D*~ D° K + decays are selected by requiring \AE\ < 25 MeV. The toes 
spectrum for these events is shown in Fig. || along with a fit with the sum of a Gaussian function 
describing the signal and an ARGUS function describing the background. The number of signal 
events is found to be 29.6±7.2. After correcting for the selection efficiencies and for the intermediate 
D° and D* + branching fractions [jl3|, the preliminary branching fraction for B° — » D*~ D° K + is 
found to be 

B{B° -> D*-D°K + ) = (2.8 ± 0.7 ± 0.5) x 10" 3 , 

where the first error quoted is statistical and the second systematic. The breakdown of the various 
contributions to the systematic error is given in Table Q]. 

Events containing B° — * D*~D*°K + decays are selected by requiring \AE + 154| < 60 MeV. 
The average position and width of AE for B° — > D*~ D*° K + is found to be in good agreement 
with expectations from B° — » D*~ D*° K + signal Monte Carlo studies. The toes spectrum of the 
selected events is shown in Fig. [/] along with a fit with the sum of a Gaussian and an ARGUS 
background function. The number of signal events found is 80.2 ± 15.3. 

To extract the B° — > D*~ D*° K + branching fraction, the contamination from decays B + — > 
D*~D* + K + , where the ii + from the D* + is not reconstructed, needs to be subtracted. This con- 
tribution has been estimated by performing the B° —> D*~D*°K + analysis on B + — > D*~ D* + K + 



signal Monte Carlo, assuming the B + — » D* D* + K + branching fraction presented in Section 5.3 



The B + —> D* D* + K + background contribution is shown in Fig. |asa small Gaussian on top of 
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the combinatorial background shape; it is estimated to be 20.6 ± 9.7 events. After subtracting this 
contribution, the preliminary B° — > D*~ D*° K + branching fraction is determined to be: 

B{B° -► D*-D*°K + ) = (6.8 ± 1.7 ± 1.7) x 10~ 3 

where the last uncertainty is systematic. The breakdown of the various contributions to the sys- 
tematic error is given in Table 111. 

5.5 B° -> D^-D^ + K° s 

The mES distribution for events reconstructed in the channels B° — » D^*' D^*' + K s is shown in 
Fig. ||. For modes involving D s, at least one decay D° — > i^7r was required. The fitted number of 
signal events is 10.1±3.7 with an estimated background of 3.4 events. The probability that the signal 
is a fluctuation of the background is 1.4 x 1CP 5 (> 5a). Most of the signal is due to the channels 
#0 _> D* + D~K S (4.7 ± 2.2 events with a background of 1 event) and B° -» D* + D*-R s (4.8 ± 2.2 
events with a background of 0.3 event). As pointed out in |l4J], the channel B° — > D* + D*~ K s is 
a CP conjugate state that could be used for sin 2/3 measurements. However, given the presently 
observed rate for reconstructing events, large improvements in the selection efficiencies are still 
needed before challenging the "golden" channels B° — » D* + D*~ as suggested in [14]. 



6 Summary 

Using about 23M BB events, we have observed several hundred completely reconstructed B 
D^*'D K decays. The following preliminary branching fractions have been measured: 

B(B° -» D*-D°K + ) = (2.8 ± 0.7 ± 0.5) x 10~ 3 

B{B° -► D*-D*°K + ) = (6.8 ± 1.7 ± 1.7) x 10" 3 

^ _ (a H 

-1.9 



in good agreement with the CLEO measurements B(B° -> D*~D°K + ) = (4.5±l[l ± 0.8) x 10 -3 



and B(B° -> D*-D*°K+) = (13.0±|;| ± 3.6) x 10" 3 ^j. 

We have observed an excess of 8.2 ± 3.5 events over a background of 1.7 events in the color- 
suppressed decay mode B + — > D* + D*~ K + , where no significant number of candidates has been 
previously seen. The corresponding preliminary branching fraction is measured to be 

B{B + -» D*-D*+K + ) = (3.4 ± 1.6 ± 0.9) x 10" 3 

Finally, several candidates have also been observed in the CP conjugate states B° — > I)(*)+Z)(*) _ Kg. 
This study confirms that the transitions b — > ccs can proceed through the decays -B — ► D^*'D K. 
To quantify more precisely this statement, we intend to measure all the individual B — > D^*'D K 
branching fractions and study the decay kinematics of these decays in the near future. 
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Figure 1: DDK decays proceed through external only diagrams (top), internal only diagrams (2nd 
line) and both (last lines) 
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Figure 2: toes an d AE 1 distributions for the sum of all neutral modes 
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Figure 3: toes an d AE distributions for the sum of all charged modes 
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Table 1: Breakdown of the various contributions to the relative systematic uncertainty on the 
B+ -► D*-D*+K + , B° -> D*-D°K + and B° -► D*-D*°K+ branching fraction measurements. 





B+ _> D*~D* + K+ 


B u _^ D*~D IJ K + 


B° -► D*-D* l> K + 


Source 


error(%) 


error (%) 


error (%) 


Tracking + Neutral efficiency 


9.7 


8.8 


8.8 


Vertexing efficiency 


10 


5.6 


8.3 


PID efficiency 


9 


5.3 


5.3 


AE requirements 


2 


7.7 


2.4 


D meson mass requirements 


13.4 


- 


- 


Intermediate BF 


5.6 


5.6 


7.5 


Background shape 


- 


4.9 


2.9 


Monte Carlo statistics 


16 


3.5 


4.3 


N £B 


1.6 


1.6 


1.6 


B+ _> D*-D*+K+ bkg 


- 


- 


19.4 


Total 


27 


16.3 


25.4 






a 
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w FS (GeV/c ) 



5.300 



Figure 4: i? + — ► D* D* + K + tties distribution 
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Figure 5: Distribution of AE versus toes f° r D*~ D° K + combinations in the data. The signal 
boxes are defined by a ±3a requirement on toes- The box \AE\ < 25MeV corresponds to B° — > 
D*~ D° K + decays, while the box \AE + 154| < 60MeV corresponds dominantly to decays B° — > 
D*~D*°K + (The 7r° or 7 from D*° — > D°7r° or D°7 is not reconstructed here) 
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Figure 6: Distribution of to E s for D* D°K + combinations with \AE\ < 25MeV. An ARGUS 
background function is used together with a Gaussian for the signal shape to fit the data. 
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Figure 7: Distribution of t«es f° r D*~D°K + combinations with \AE + 154| < 60MeV (B — > 
D*~ D*° K + signal region). An ARGUS background function is used together with a Gaussian for 
the signal shape to fit the data. The B + — > D*~ D* + K + background contribution is shown as a 
small Gaussian on top of the combinatorial background shape. 
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Figure 8: B° -» D^'-DW+Kg m E s distribution 
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